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Abstract 

We  investigate  the  compression  of  rubber  blocks  bonded  to  rigid  metal  plates  for 
several  different  geometries,  namely  for  rectangular  and  V-shaped  blocks.  We  obtain 
reasonably  simple  solutions  in  series  form  relating  the  force  and  the  displacement  using 
the  ’’pressure  method”  developed  by  Gent.  Our  initial  assumption  concerning  the 
incompressibility  of  the  rubber  is  subsequently  relaxed  to  ’’near  incompressibility”  to 
yield  more  accurate  solutions.  We  demonstrate  a  very  good  agreement  between  the 
theoretical  solutions  and  the  results  of  linear  finite  element  simulations. 

Key  words:  elasticity,  compression,  incompressibility,  near  incompressibility,  partial  differ¬ 
ential  equations,  finite  element  simulations 

1  Introduction 

In  this  note  we  examine  the  compression  of  rubber  blocks  bonded  to  rigid  metal  plates  for 
several  specific  geometries.  Rigorous  analysis  of  such  components  is  nontrivial  since  the 
solutions  are  complicated.  However,  for  design  purposes  it  is  most  desirable  to  produce 
(if  possible)  reasonably  simple  closed  form  solutions  for  these  geometries  that  provide  a 
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good  approximation  to  the  actual  behavior  of  the  components.  A  number  of  approximate 
methods  have  been  developed  using  simplifying  assumptions  such  as  incompressibility  of  the 
rubber  and  linear  elastic  behavior.  Here  we  adapt  an  approach  developed  by  Gent,  [2,  3,  4], 
called  the  ’’pressure  method”.  When  a  bonded  block  is  compressed,  the  deformation  can 
be  considered  to  take  place  in  two  stages:  a  simple  homogeneous  compression,  (F\),  which 
is  equivalent  to  compression  between  lubricated  plates  and  subsequent  shear  deformations, 
(F2),  that  restore  points  in  the  planes  of  the  bonded  surfaces  to  their  original  positions. 
Hence,  assuming  superposition,  the  compression  force  is  given  by  F  —  F1  +  F2.  Solutions 
derived  from  the  pressure  method  approach  have  been  quite  successful  and  are  often  quoted 
in  rubber  engineering  design  handbooks  ([6]).  As  we  described  in  [1],  solutions  using  this 
method  give  very  good  approximations  to  the  compressive  force  in  the  case  of  rectangular  and 
annular  blocks.  As  we  shall  demonstrate,  the  assumption  of  incompressibility  can  be  relaxed 
and  a  good  agreement  can  be  established  between  the  resulting  theoretical  formulation  and 
the  numerical  finite  element  results  for  slightly  compressible  rectangular  and  annular  blocks. 

In  this  paper  after  a  general  introduction  to  the  pressure  method  we  revisit  the  case  of  a 
rectangular  bonded  block  and  present  in  detail  the  results  for  the  slightly  compressible  case. 
Next  we  apply  the  pressure  method  to  a  V-shaped  mount  and  extend  the  solutions  to  the 
slightly  compressible  case.  At  the  end  of  this  note  we  present  a  comparison  of  theoretical 
and  finite  element  analysis  predictions. 


2  The  pressure  method 

In  this  section  we  provide  a  brief  introduction  to  the  pressure  method  and  present  the  general 
strategy  used  in  obtaining  the  governing  equations  for  different  geometries.  Other  descrip¬ 
tions  of  the  pressure  method  can  be  found  in  [2,  3,  4],  As  stated  in  the  introduction  we  use 
the  assumption  that  the  shear  modulus  is  constant  (i.e.,  there  is  a  linear  relationship  between 
stress  and  strain)  and  thus  our  results  apply  for  small  strains.  For  small  compressions  the 
load-deflection  relation  may  be  written  as: 

J  =  Ef(S)s, 

where  A  is  the  cross-sectional  area,  f(S)  is  the  shape  function,  S  =  J^ondedaxea  i  E  is  Young’s 
modulus  for  the  rubber  and  e  is  the  average  compressive  strain.  The  goal  from  the  designer’s 
point  of  view  is  to  compute  f(S)  for  different  geometries.  The  deformation  is  considered  to 
occur  in  two  stages:  first  there  is  a  simple  homogeneous  compression  and  then  subsequent 
shear  deformations  restore  points  in  the  planes  of  the  bonded  surfaces  to  their  original  posi¬ 
tions.  We  denote  the  compressive  force  resulting  from  the  simple  homogeneous  compression 
by  F1.  The  compressive  stress,  cq,  for  the  first  deformation  is  uniform  and  given  by: 
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when  displacements  occur  equally  in  the  x,  y  plane,  or 


Fi  Ee 

A  1  1  —  v2  ’ 


where  v  is  Poisson’s  ratio,  when  the  displacements  are  restricted  to  one  direction  (e.g.,  long 
strip).  For  incompressible  materials,  v  is  0.5  so  that  =  §•  Then  for  ^  =  Efi(S)e, 

fi(S)  ranges  from  1  to  |. 

An  empirical  form  for  the  fi(S)  term  (see  [3])  is: 


fi(S) 


4  2  ab  + 12 

3  ~~  3  a?  +  b2  +  2 12  ’ 


(2.1) 


where  2a,  2b  are  sides  of  the  rectangular  cross-section  or  axes  of  the  elliptic  cross-section. 
Note  that  this  formula  gives  1  when  displacements  occur  equally  in  the  x,y  plane  (i.e. ,  square 
block)  and  equals  |  for  an  infinite  block. 


p+  4k  dx 
dx 


Figure  1:  Forces  on  an  element 


According  to  Gent,  under  certain  conditions,  the  stresses  required  to  maintain  the  second 
deformation  can  be  deduced  readily  by  analogy  with  Reynolds’  theory  of  thin  lubricating 
films  of  a  viscous  liquid.  For  relatively  thin  blocks  the  stresses  are  represented  in  the  form 
of  a  hydrostatic  pressure  P  that  is  maximal  in  the  center  of  the  block  and  falls  to  zero  at 
the  free  surfaces.  The  governing  relation  for  P  can  be  derived  for  different  geometries.  To 
illustrate  the  method  we  present  the  derivation  of  P  in  detail  for  a  thin,  bonded  infinite 
strip. 
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Consider  static  equilibrium  of  forces  on  the  element  in  Figure  1  where  P  denotes  the  pressure 
and  a  denotes  the  shear  stress  in  the  x,  y  plane: 


dxdydz  =  — — dzdxdy , 
ax  dz 


(2.2) 


and  thus 


dP  da 
dx  dz 


(2.3) 


Since  the  pressure  is  assumed  constant  through  the  thickness  ^  is  independent  of  z,  which 
implies  that 


d2a 

dz2 


=  0. 


(2.4) 


The  shear  stress,  a,  is  given  by 


a  =  -G 


dk 
dz  ’ 


(2.5) 


where  k  is  the  x-displacement  of  the  volume  element  and  G  is  the  shear  modulus  of  the 
rubber 

From  (2.4)  it  follows  that 


d3k 
dz 3 


=  0 


and 


dP  _  d2k , 
~di  ~ 


(2.6) 


(2.7) 


Equation  (2.6)  indicates  that  an  initially  plane  cross  section  takes  a  parabolic  shape  in  the 
deformed  state,  with  a  profile  given  by: 


km  Vi/’ 


(2.8) 


where  k  =  0  when  z  =  ±|  and  k  —  km  when  z  =  0.  Substituting  (2.8)  into  (2.7)  we  obtain: 


dP  _ 


(2.9) 


To  find  km  we  equate  the  volume  of  the  bulge  with  the  volume  displaced  by  the  compres¬ 
sion  of  the  block  lying  within  the  plane  at  x  (see  Figure  2).  Only  a  unit  length  is  considered 
in  the  case  of  an  infinite  block. 
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Figure  2:  Volume  of  the  bulge 


The  change  in  volume  is  given  by 

AV 


On  the  other  hand, 


etAx. 


AV 

~V~ 


P 

Kr 


where  K  is  the  bulk  modulus  of  the  material.  Thus 


AV 

~v 


2  Akmt 


etAx 


3  Axt  Axt 


P 

K' 


which  yields 


2  A  kr 


3  Ax 
dkm 


dx 


—  e  — 


K 


h  _  p \ 

2  £  K^' 


(2.10) 


Note  that  if  K  —  oo  (incompressible  material)  then  | e .  Combining  (2.9)  and  (2.10) 

we  obtain 


d2P  G .  P, 

-  ~12fi(e  -  K>- 


(2.11) 


When  deformation  occurs  in  the  y-direction  as  well  we  have 


d2P  d2P  _  G  _P_ 
dx2  dy2  t 2  K 


(2.12) 
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In  objects  with  cylindrical  symmetry  we  obtain 


cPP  1  dP  _  G  P 

dr2  ^  r  dr  t 2  K 


(2.13) 


If  incompressibility  is  assumed  then  the  equation  is  the  same  with  the  term  E  removed  (i.e., 
take  K  —  oo).  Once  the  pressure  P  is  found,  the  corresponding  compressive  force  can  be 
obtained  by  integration  and  /2(5)  (the  component  of  the  shape  function  arising  from  the 
second  step  of  the  deformation)  can  be  determined.  Thus  for  a  long  bonded  strip  the  solution 
of  (2.11)  with  boundary  conditions  P  —  0  at  x  —  ±w/2,  and  dP/dx  —  0  at  x  —  0  yields 


P  = 


6Ge 

~P~ 


in  the  incompressible  case.  Thus 


F2  = 


Pdx 


Gew 3 
t2 


(2.14) 


(2.15) 


For  the  long  bonded  strip  the  force  per  unit  length  associated  with  the  simple,  homogeneous 
compression  is  given  by 

Fi  =  \Eew,  (2.16) 

where  E  is  the  Young’s  modulus  of  the  rubber.  Thus, 

4  (  w2  \ 

F  =  Fi  +  F2  =  -Eew  ( 1  +  —  j  ,  (2.17) 


and  since  the  shape  factor  for  the  long  bonded  strip  is  S  —  |£,  we  obtain  that  the  shape 
function  is 

/(S)  =  |(l  +  S2). 

As  we  noted  in  [1]  Equation  (2.12)  is  analogous  to  the  equation  governing  the  torsion  of 
a  prism.  Solutions  of  the  torsion  problem  exist  for  many  different  geometries  [7,  8,  9]  and 
can  be  used  directly. 

Historically,  the  general  incompressible  problem  was  solved  by  Gent  and  Meinecke  in 
1970  [3].  In  1994  Gent  extended  his  solutions  for  the  infinite  strip  (1-D)  and  the  disk 
(axis-symmetric  geometry)  to  a  compressible  material  [4],  The  axis-symmetric  solution  was 
generalized  by  Kelly  to  the  annulus  [5].  Our  work  extends  the  incompressible  solution  for 
rectangular  blocks  to  compressible  material  and  uses  this  result  to  find  solutions  for  V- 
shaped  rubber  blocks.  We  also  compare  the  theoretical  results  with  numerical  simulations 
using  finite  element  methods. 
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3  Rectangular  blocks 


In  this  section  we  review  the  compression  of  bonded  rectangular  blocks  of  length  /,  width 
w  and  thickness  t.  In  [1]  we  outlined  the  results  for  the  incompressible  and  the  slightly 
compressible  case;  here  we  present  the  analysis  in  detail  for  the  latter  case.  Results  for 
the  former  can  be  found  in  [1,  3],  but  we  note  that  they  are  not  expected  to  provide  a 
good  approximation  unless  some  compressibility  correction  is  used  especially  when  the  shape 
factor,  S,  is  large. 

The  compressive  force  F1  corresponding  to  the  pure  homogeneous  deformation  is  given 
by 

§  =  Eeh(S), 

where  fi(S)  is  the  empirical  factor  given  by  (2.1),  with  a  =  1/2,  and  b  =  w/2.  The  cross- 
sectional  area  A  is  Iw.  As  noted  in  [1]  for  blocks  of  moderate  to  large  shape  factors  F1  is 
much  smaller  than  F2,  so  uncertainty  in  fi(S)  is  usually  not  critical. 

We  derive  the  force  F2  by  analogy  to  Reynolds’  squeeze  film  theory.  In  the  slightly 
compressible  case  the  pressure  satisfies  the  equation 


d2P  d2P  _  12  Gf  P\ 

dx2  dy2  t2  v  Kj  ’ 

with  boundary  conditions 


P  =  0  at  x  =  ±w/2,  y  =  ±£/ 2,  and 

dP  n  dP  n 
—  =  0,  — -  =  0  at  x  =  y  =  0, 
ox  oy 


(3.18) 


(3.19) 

(3.20) 


where  K  is  the  bulk  modulus  of  the  rubber.  Since  the  pressure  is  symmetric  with  respect  to 
the  y-axis  and  P  =  0  at  x  =  ±w/2  we  propose  to  find  a  solution  to  (3.18)  in  the  form 


P{x,y) 


ni ra: 
an  cos  —  Yn, 


n=  1,3,5,... 


W 


(3.21) 


where  an,  n  =  1,3,5,...  are  constant  coefficients  and  Yi,Y3,...  are  functions  of  y  only. 
Substituting  (3.21)  into  (3.18),  formally  we  obtain 


which  yields 


00  n27T2 


E 


n=  1,3,5,... 

12G 


o 

W2 


TITT  X  r 

COS  In  /  4 

n= 1,3,5,... 


nnx 

a„  cos - Y„ 


w 


w 


t2 


£  — 


K 


E 


mrx , 


(Ir,  cos  ■ 


-Yn 


n= 1,3,5,... 


W 


E 

n= 1,3,5,... 


U'KX  , 

an  cos - Yn 

w 


n27T2 

^ n  2 
W 2 


U'KX 


■  COS 


+ 


w 


12  G 
¥k 


CLr j  COS 


U'KX  \ 

I  in 
W  J 


(3.22) 


12  Ge 
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We  expand  the  right  side  in  Fourier  series  for  —w/2  <  x  <  w/2  to  obtain 


12  Ge 


JS,  12 Ge  4  ,  n— i  rara: 

—  >  — ^ - (  — 1)  2  cos 

^  t2  m r v  ' 


n= 1,3,5,... 

Thus  Yn(y)  satisfies  the  following  ordinary  differential  equation 


w 


Y" 

n 


'n2  7T2  12G\ 

w2  t2K ) 


Yn  =  -Uf  4  (-1)^, 

tz  mra-T, 


with 


Yn(±£/2)  =  0  and 

^(o)  =  o. 


The  general  solution  of  (3.24)  is 


^  ,  ,  .  .  .  .  .  1  12 Ge  4  .  .  n^i 

Yn(y)  —  C\  smh  Xny  +  c2  cosh  Xny  +  — — - (-1)  2  , 


A2  t 2  mran 


(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 


with 


A 


n27T2 


+ 


12  G 
¥k' 


Since  the  pressure  is  symmetric  with  respect  to  the  x-axis,  it  follows  that  c\  =  0.  From  (3.25) 
we  obtain  that 

1  12 Ge  4  ,  1 


c2  = 


-(-I)' 


A2  t2  mran  v  '  cosh  AA 


We  note  that  in  this  case  (3.26)  is  also  satisfied.  Thus 


12  Ge 


E 


with 


t2 


Ad  = 


(-1) 


n—1 

2 


1  - 


cosh(Anr/) 

cosh(Anf/2) 


mrx 


cos 


w 


n27T2 


+ 


12  G 
¥k' 


We  integrate  P  over  the  metal  plate  to  obtain 

r*/ 2  r»/ 2  ,  24 Ge 

Pdxdy  = 

I -i/2  J-w/2 


rt/2,  rw/2 

F2  =  /  /  Pdxdy  = 

J-e/2  J-w/2 


E 


4w£ 


8w  nh  Kl 

f2  n=ts,5,...  n2^2X2n  n%2A3  &n  2 


This  series  converges  rapidly  and  can  be  computed  easily  for  different  ratios  of  w/l. 
Adding  F\  and  F2  yields 


(3.28) 


(3.29) 


F 

A 


F\  +  F2  _  1 
A  “  A 


24Ge  “  4  w£  8w  XJ 

t2  n%2An  n2n2\l  an  2 


.  (3.30) 
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This  result  was  compared  to  numerical  simulations  using  linear  finite  element  methods  in 
[1],  The  half-width  of  the  block  (w/2)  was  kept  constant  at  0.0254  m  (1.0  in)  while  the  half 
length  {lj 2)  and  the  thickness  t  were  varied.  The  material  was  assumed  to  be  linearly  elastic 
with  Young’s  modulus  of  4.137  •  106  Pa,  (600  psi)  and  a  Poisson’s  ratio  of  0.4995  (slightly 
compressible  material).  The  theoretically  and  numerically  predicted  vertical  force  needed  to 
compress  the  rectangular  blocks  are  depicted  in  Figure  3.  Discrepancy  between  the  finite 
element  results  and  theory  is  less  than  10%  which  is  very  satisfying  considering  that  both 
methods  are  approximate. 
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Comparison  of  FEA  results  for  F  with  theory  for  rectangular  blocks 
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Figure  3:  Comparison  of  FEA  results  for  F  with  theory  for  rectangular  blocks 


4  V-shaped  blocks 

In  this  section  we  derive  the  total  axial  force  corresponding  to  a  small  vertical  displacement  5 
for  V-shaped  blocks  under  different  assumptions.  V-shaped  rubber  blocks  are  used  in  rubber 
springs  to  achieve  particular  combinations  of  stiffnesses  in  the  three  orthogonal  directions. 
They  can  be  found  for  example  in  railway  bogie  suspensions.  Our  goal  is  to  develop  expres¬ 
sions  for  the  shape  function  to  aid  in  the  design  of  such  blocks  by  predicting  the  stiffness 
parameters. 
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4.1  Infinitely  long  V-shaped  blocks 

We  start  by  analyzing  the  case  of  infinitely  long  Chevron-shaped  blocks.  We  assume  that 
the  length  of  the  block  is  much  larger  than  its  width,  thus  we  can  neglect  the  displacements 
at  the  two  ends.  We  also  assume  in  this  derivation  that  the  rubber  is  incompressible.  The 
cross-section  of  the  V-shaped  block  is  depicted  in  Figure  4.  Since  the  block  is  symmetric, 
we  consider  only  one  leg  of  the  Chevron-shaped  mount  first. 


Figure  4:  Cross-section  of  the  infinitely  long  V-shaped  block 


The  vertical  displacement,  <5,  may  be  resolved  into  its  components  Ssin9,  which  is  normal 
to  the  metal  plate  and  causes  compression  in  the  rubber,  and  ScosO,  which  is  tangential  to 
the  metal  plate  and  causes  shear.  We  want  to  calculate  the  vertical  force  per  unit  length, 
Fa,  on  one  leg  of  the  chevron.  This  is  given  by 

F  =  F  -t-  F 

J  a  c  '  x  si 

where  Fc  and  Fs  are  the  contributions  (per  unit  length)  from  the  compression  and  shear 
deformations  respectively.  First  we  consider  compression,  which  is  assumed  to  take  place 
in  two  stages  as  explained  in  the  introduction.  The  force,  F1,  associated  with  the  simple, 
homogeneous  compression  is  in  the  form  given  by 

Fi  _  4  6  sin  9  _  4 G5  sin  9 
~A  ~  3  t  “  t  ’ 

with  the  cross-sectional  area  A  =  w/2.  Since  incompressibility  is  assumed,  we  have  E  —  3 G. 
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To  compute  F2  we  again  use  the  analogy  to  Reynolds’  squeeze  film  theory.  As  we  showed 
in  Section  2  the  local  hydrostatic  pressure  P  sustaining  the  bulge  profiles  solves  the  ordinary 
differential  equation 


where  km  is  the  maximal  bulge  (see  Figure  2).  The  volume  of  the  bulge  is  approximated 
by  V  ~  \kmt}  while  the  displaced  volume  per  unit  length  is  given  by  V  —  xS  sin  9.  Hence 
km  =  ^xS  sin  9. 


2 t 

Thus, 


and  the  solution  is 


dP  12  GxS  sin  9 

dx  t 3 

„  w 

P  —  0  at  x  =  ±— 

2 


6G<5  sin  9  ( w2 


The  pressure  acts  on  the  metal  plate  and  therefore  the  corresponding  force  per  unit  length, 
F2,  is  obtained  by  integrating  P  along  the  plate 


Fo  = 


" /2  GGS  sin  9  ( w2 


—  x2  \  dx  — 


G5  sin  9wz 


This  implies  that  the  compression  force  acting  perpendicular  to  the  metal  plate  is  given  by 

_  _  _  4  5  sin  9w  G6  sin  9w3 

F  =  F1  +  F2  =  -E - 1 - - - 

1  2  3  2 1  2t3 

4  5 sin 9  (  w2\ 

=  PE—{l  +  w)’ 

and  the  axial  component  Fc  of  the  compression  force  F  is 

4  8  sin2  9  (  w2\ 

c  3  t  \  At2) 

To  find  the  total  axial  force  Fa  we  also  need  to  compute  the  vertical  component  of  the  shear 
force  induced  by  the  tangential  component  of  the  displacement.  The  shear  force  can  be 
determined  as 

GAS  cos  9 
F,  =  — t — 1 

and  its  contribution  to  the  axial  force  is 


Fs  =  Ft  cos  9  = 


GAS  cos2  9 
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Thus  the  axial  force  on  one  leg  of  the  V-shaped  block  is 


GA5  cos2  9  4  <5  sin2  9  f  w2\ 

Fa  =  F,+Fc=  -  +3^— ^(l+4id' 


which  implies  that  the  total  axial  force  per  unit  length  is 

Fa  =  2  Fa  = 


Gw5  cos2  9  4  „  5  sin2  9  /  w2 


t 


H — Ew- 
3  t 


1  +  4F  ’ 


since  A  =  w/2.  It  is  readily  seen  that  this  expression  reduces  to  our  result  for  the  infinite 
strip  when  8  =  7t/2. 


4.2  Finite  V-shaped  blocks 


Figure  5:  Finite  V-shaped  block 


As  we  did  in  the  previous  example,  we  consider  only  one  leg  of  the  finite  V-shaped 
block  of  length  l.  As  before,  the  vertical  displacement  5  induces  compression  and  shear  in 
the  rubber  block.  We  first  find  the  axial  component  of  the  compression  force.  Recall  that 
F  =  Fi  +  F2,  where  Fi  comes  from  the  homogeneous  compression.  We  use  the  empirical 
factor  (2.1)  proposed  by  Gent  and  Meinecke  [3]  to  obtain 


F\ 


ES  sineA  (i 

2 1  13 


2  w£/4.  +  t2  \ 

3w2/4  +  £2/4  +  2 12 )  ' 
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F2  is  again  derived  by  Reynolds’  theory.  First  we  have  that  the  pressure  P  satisfies  the 
following  equation: 


d2P  d2P  G  _  .  . 

dx2  '  dy2  12tzdsm  ’ 

(4.31) 

7-,  W  l 

P  =  0  at  x  =  — ,  y  =  ± 

(4.32) 

dP 

——  =  0  at  x  —  0. 
dx 

(4.33) 

Let  us  denote  the  right  side  of  (4.31)  by  M.  Then  this  equation  is  analogous  to  the 
torsion  problem  described  in  [7],  and  very  similar  to  the  compression  of  a  finite  rectangular 
block.  Thus  we  can  use  arguments  similar  to  those  in  Section  3  to  obtain  a  solution.  The 
boundary  conditions  at  x  —  and  x  —  0  are  satisfied  by  looking  for  P  in  the  form  of  a  series 


p=  E 


mrx 

cos - Yn, 

w 


n= 1,3,5,... 

where  Yi,Y3,  . . .  are  functions  of  y  only.  We  can  argue  as  in  Section  3  to  find  Yn  and  P. 
With  routine  calculations  we  thus  obtain 


P  ~  8 sin  9  f  -±^(-1) 


f3 


777T 

1,3,5,...  nn  An 


n—  1 
2 


1  - 


cosh (A ny) 
cosh(An^/2) 


cos 


m^x 

j 

W 


where  A 1  = 


n  w2  ■  To  find  P2  we  integrate  P  along  the  bonded  plate  of  length  i 

1/2  rw/2 


/1/2  fW/2 

/  Pdxdy 
-a  2  Jo 


>-1/ 2 

12G 

~p~ 


S  sin  9  Y. 


n=  1,3,5,... 


Awl  8  W  A  nl 

-  n%2A3  tanh  ~2~ 


(4.34) 


F2  = 


2ts 


G8  sin  9wzl  I  1  — 


192  w 


7 r 


5  l 


>oo 

71=1,3,5 

J_  _ 
,•••  n4 

7T4 

—  we 

96 

oo 

1 

ukI  \ 

E 

1,3,5,... 

—  tanh 

n5 

2w  J 

This  series  converges  rapidly  and  it  can  be  computed  easily  for  different  ratios  j.  For 
example  if  w  —  l  (which  is  not  a  square  cross  section  of  one  leg  in  our  notation)  we  obtain 

,^,sm,4 

For  w/2  —  l  (which  is  a  square  cross  section  of  one  leg  in  our  notation)  we  find 

j-,  0.1716  .  y 

F2  =  ~Kfi~GSs,n9Y' 
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In  the  general  case  the  compression  force  is  given  by 


F  =  F1  +  F2  =  Ed  sin  9 


w£  (4 


w£/A  + f2 


2 1  \3  3  w2/A  +  £2 /A  +  2f2 


+ 


2f3 


G5  sin  9ws£  1  — 


192  w 


7 r 


5  £ 


E 


n=l,3,5,... 


1  mv£ 

— tanh  - 

n5  2  w 


E-—6  sin  6  1  —  —  — 
2 1  "  ‘ 


w£/ A  +  f2 


+ 


3  3  «;2/4  +  £?  j  A  +  2f2 

192  w 


7 T5  l 


^  1  7l7lT 

2^  -5 tanh 


n=l,3,5,. 


rr 


2tc 


The  axial  component  of  the  compression  force  is  Fc  =  F  sin  9. 

The  shear  force  induced  by  the  tangential  component  of  the  displacement  is  given  by 


Ft  = 


GwS  cos  9 
2 1 


and  its  axial  component  is  Fs  =  Ft  cos  9.  Thus  we  can  evaluate  Fa  =  Fc  +  Fs  for  different 
sized  blocks.  Again  the  total  axial  force,  Fa  =  2 Fa,  reduces  to  our  result  for  rectangular 
blocks  when  9  =  tt/ 2. 


4.3  Extension  to  slightly  compressible  V-shaped  blocks 

In  the  above  analysis  we  made  the  assumption  that  the  rubber  is  incompressible.  In  reality 
rubber  is  slightly  compressible,  with  a  large,  but  finite  bulk  modulus.  The  usual  incompress¬ 
ibility  assumption  stems  from  the  fact  that  the  shear  modulus  of  the  rubber  is  relatively  low, 
about  a  thousand  times  smaller  than  its  bulk  modulus.  Thus  rubber  accommodates  stresses 
in  a  way  such  that  the  change  in  volume  is  negligible  compared  to  its  change  in  shape. 

In  this  case  we  need  to  modify  our  results  for  the  F2  component  of  the  compression  force. 
The  other  forces  are  not  affected  by  the  bulk  modulus.  The  equation  we  must  solve  in  this 
case  is  the  following 


d2P  d2P 
dx2  dy2 


-Upsine+^P, 


w  £ 

P  =  0  at  x  =  — ,  y  =  ±-, 


dP 

dx 


0  at  x  =  0, 


(4.35) 

(4.36) 

(4.37) 


where  K  is  the  bulk  modulus  of  the  rubber.  Let  us  again  denote  the  first  term  on  the  right 
side  of  (4.35)  by  M.  We  can  solve  the  equation  in  the  same  manner  as  we  solved  (4.31).  The 
boundary  conditions  at  x  =  ^  and  x  =  0  are  satisfied  by  seeking  P  in  the  form  of  a  series: 


p=  E 

«=1,3,5,. 


mvx 

cos - Yn, 

w 
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where  Yi,  Y3, . . .  are  functions  of  y  only.  Arguing  as  before  we  obtain  that 


f3  „  ft*  nn X2nK  ’ 

n— 1,0,0,...  n 


n—  1 
2 


1  - 


cosh  A ny 
cosh  An| 


mrx 


cos 


w 


with 


a: 


n2  7T2  12G 

+ 


n  9 


t2K 


To  find  F2  we  integrate  P  along  the  bonded  plate 

rt/ 2  rw/2 


/l/Z  f'W/Z 

/  Pdxdy 

-1/2  JO 


1-1/ 2 

12  G 

I5” 


<5  sin  9  y. 


4  wl 


8 w  .  ,  A„£ 


n=l,3,5,... 


n2n2X2  n2TT2Xi  2 


Once  the  force  F2  is  obtained  the  remaining  analysis  is  the  same  as  in  the  previous  subsection 
since  the  other  forces  are  not  affected  by  the  bulk  modulus  and  the  assumption  that  the 
rubber  is  slightly  compressible. 


4.4  Finite  element  analysis  results 

Linear  finite  element  analysis  was  performed  using  ANSYS  5.6  to  verify  the  theoretical 
solution  for  the  finite  V-shaped  block.  3D  analysis  was  performed  using  Solid95  elements 
with  reduced  integration.  The  material  was  assumed  to  be  linearly  elastic  with  Young’s 
modulus  of  4.137  •  106  Pa,  (600  psi)  and  a  Poisson’s  ratio  of  0.4995  (slightly  compressible 
material).  To  cover  a  reasonable  range  of  shapes,  the  angle,  9  was  fixed  at  7t/4  and  7t/3. 
The  width  w  was  kept  at  0.0505  m,  (2  in),  while  the  thickness  varied  from  0.0031  to  0.0126 
m,  (0.125  to  0.5  in),  and  the  length  varied  from  0.0126  to  0.203  m,  (0.5  to  8  inches).  A 
y-displacement  of  2.53990  •  10-5  m,  (0.001  in)  was  imposed  on  the  top  surface  to  simulate  a 
small  compression  S.  Figures  6-7  show  the  comparison  of  the  finite  element  and  the  theoretical 
results  for  the  vertical  force  needed  to  compress  V-shaped  blocks  of  different  angles  and  sizes 
by  5  in  the  slightly  compressible  case.  There  is  a  very  good  agreement  between  the  two 
results  with  a  discrepancy  of  less  than  10  %.  This  is  quite  satisfactory  considering  that  both 
methods  are  approximate.  The  largest  discrepancy  occured  when  the  length  i  was  small  and 
the  width  t  was  large.  This  is  not  surprising  since  in  that  case  some  of  the  assumptions  of 
the  pressure  method  do  not  apply. 


5  Acknowledgments 

We  are  grateful  to  Gail  G.  Licayan  for  her  contributions  in  organizing  some  of  the  material  for 
this  paper.  This  research  was  supported  in  part  (HTB  and  GAP)  by  the  Air  Force  Office  of 
Scientific  Research  under  grants  AFOSR  F49620-98-1-0180  and  AFOSR  F49620-01-1-0026. 


15 


Comparison  of  FEA  results  for  F  with  theory  for  0=ji/4 
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Figure  6:  Comparison  of  FEA  results  for  Fa  with  theory  for  6  =  ir/ 4 
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Figure  7:  Comparison  of  FEA  results  for  Fa  with  theory  for  9  =  7t/3 
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